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ABSTRACT 


This paper presents an independent verification of the 
«location method as a technique for calculating the lift 
on an oscillating airfoil in an unstaggered cascade 
immersed. in transonic flow. This method was originally 
proposed by Gorelov. Results presented here differ some- 
what from those presented by him. Two formulations are 
shown; one is purely numerical, the second employs an 


analytic expansion for small frequency. 
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I: INTRODUCTION 


The analysis of unsteady transonic flows in aircraft 
turbopropulsion is an area of intense current interest. 
Rising fuel prices and increasing thrust requirements both 
point toward the need of turbomachinery capable of performing 
well with transonic or supersonic internal flow. But, 
increased flow has increased both the costs and uncertainties 
of engine designs.  Flutter problems have already become a 
major consideration in engine development. Problems unfore- 
seen in earlier days of turbine engine production have 
caused long development delays, or forced acceptance of 
engines producing less than their initial design thrust. 
These uncertainties cannot be avoided when an attempt is 
made to extend the state of the art, but they can be reduced 
by extending the range of analytical modeling. 

Such extension must now be done piecemeal. The three- 
dimensional flows in turbomachinery, including the simul- 
taneous effects of boundary layers, rotation, finite blade 
thickness, spanwise Mach distributions, and shocks, are 
well beyond present capability. Perhaps one day complete 
E Ss will be practical, but it is not today. The best 
that can be done now is to approach the problem from one 
aspect at a time. Flow through a two dimensional cascade 


has been a useful tool in this partial analysis. 
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This thesis was originally to have been an extension of 
the work of Elder [1] and Schlein [2] to the case of a 
Staggered cascade. Their work, based on Teipel's [3] 
linearization of the unsteady transonic small perturbation 
equation, analyzed transonic flow through oscillating 
unstaggered cascades by use of the collocation method. 
While the problem was easy to state, it was difficult to 
solve. Both Elder and Schlein had encountered difficulty 
in employing the collocation method. Therefore, it was 
Gecided that verification of the basic collocation solution 
presented by Gorelov [4] using a different linearization 
would-be a worthwhile goal in itself. 

The following investigation presents a verification 
of the development in [4], along with numerical results and 


Suggestions for further work. 
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II.  UNSTEADY TRANSONIC FLOW THEORY 


BEcnssdering inviscid Elow only, the following four 


equations govern the aerodynamic flow problem at hand: 


The equation of state 


DS pk 


and the equations for the conservation of 


l. Mass: div(oV) + 22 = 0 
2. Momentum: + : Mp. -— X9 | 
DIN Energy: == = 
where ^ 
v = velocity 
p = pressure 
S = entropy 
R = universal gas constant 
T = temperature 
t = time 
0 = density 
De ok a 
Dt üt 3x Ot day dt 


πι ο 


il =z ) 


(Eis) 


(II-4) 


The analysis starts with-a uniform flow From edm ty. 


This flow has velocity U parallel co the x- axis. This 
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formulation can be simplified by working with the total 


metocity potential, 6, where 


e 
I 
Qo 

ls 
ll 

Ον 

|| 
Qo 
νά 


USD 


X component of velocity 


Q 
νά 
νά 


Q2 
9 
QJ 
M 


v = = 0 = y component of velocity = — (II-6) 


Qo 
ke 
Ke 


Thus, the initial uniform flow is represented by the uniform 
flow potential 


ᾧ τοι ο ος (II-7) 


This notation may be applied to the conservation equations 
for mass and momentum. The equation for conservation of 


mass 


> garen $ 
αι ν py) + TE 0 (II-8) 


becomes for two-dimensional unsteady flow 








ð (pu) 9 (pv) ος 
ox oy at : 
(22 * u 28 « o t ροζ. € $T) = 0 ΕΕ 9) 
ot ox 
Iu t 
30 3p 30 . Do 
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and 


u = Ὃν and y = um 
Thus 
Do 
— + = 
a O 0 
and 
m 
+ = — e ae xk 
Pex "yy o Dt eri ο) 


Thus 


Peetying this to equation (II-10) yields 


> TOD = = =, =P (II-1l1la) 
XX YY pa E 
$ ae e = - EE Πρ. ο. υπ το) 
Rt yy s x y E 
1 
= -= Fo do) 5 4 Pel (II-11c) 
pa 
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where Vis the gradient operator, P - 2E P = s P, - 2R s 


AA E E 
Laying this aside for the moment, consider the momentum 


equation (11-3) 


Dv. il 
"Τ᾽ 
> 
DURS oy 
Dt EU (v V) V 
Thus 
UA NUI Mays (II-12) 
ot 
y = Vo 
SO 
av a 2 = ον (11-13) 
ES 5c (79) cr 
and 
2 
Vv- vV) v = VERE ccv x (ve v) 
where 
ο. ο JE 
> > 
= vV Y 
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Bor irrotational flow 


V X ν = 0 
thus 
2 
> VV 
(v * V)v = — Ga 
Thus 
+ vo, + v = 0 5 


which after integration along a streamline becomes 


2 


dp pn 
ca ο E) (11-16) 
τ. τπιεοτπ flow from infinity F(t) = 5 ο and thus 


the final result is 


2 
dp Ve A 
[= + στ. 


; es πι RUN 


bo] Fe 
E 


Differentiation with respect to t gives 


20. 





= l ον 
Ehe m pud 1 5 ο (11-18) 
Eom (II-3) 
Dv 
=Vp = 9) δε πι ο) 


Substitute (11-18) and (II-19) into (II-1llc) to obtain 


τα 2 
i l Dv I 3V 
=F = — ER en nt = 
"xx "yy ο. Dit i Per j ο τ; 
This may be further simplified 
DV 9v > -> 
pc AN UE 
ZEN στ. 
D ro. 
Hence 
φ + 0 = Lio Y lov?) + + = av’, (II-21) 
ER yy 44 dt 2 ce 2 ος 
Expanding terms 
ὃν. =O! MU. E 
(VO) (FF) = Vo (Evo) = Sep zh Br yt 
gy? $ o $ = 2 2 o QD i 
EMEND e mu cp. poe. cxx 
o 24 ® 25 29,9 qa 
= x = n ZH T Y y 


2 l 





Q> 
< 


EE 
2 3gl(V9) + (V9)] 


ho| t— 
a 
"| 

| 


X XE XE X y yt yt y 


The final result obtained by combining terms is 


2 2 


D ® <9 0 Py 
Bee? xx * (2 > TF) yy 7 2 
a a a 
dn. 29s 1 
EN 0 EE i = yp = 0 (15-22. 
κ xt a? ye a? tt 


BN-result is valid for irrotational, inviscid, two- 
dimensional, unsteady, compressible flows where gravity has 


been neglected. 
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III. SMALL PERTURBATION THEORY OF TRANSONIC FLOW 


A. GENERAL CASE 

A thin body at a small angle of attack will cause only 
a slight disturbance in the fluid. A flat plate 15 8η 
exanple.  Consicer flow past a flat plate at angle of 


aam ack a. 


(1) (2) 


The flow at (2) must be parallel to the plate. To achieve 
this, small disturbance velocities u' and v' must be added 


to the free stream velocity yielding 


The potential of the disturbed flow may be considered 
as the sum of the uniform flow potential, es = Ο.Χ, and the 


disturbance potential, 6 


6 = Ö + ὁ (III-1) 
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Thus 


9 zu are: (III-2a) 


$ = V (ITI=2D) 


Esa function of time, then 


This result may be substituted into (II-22) leading to 


- 


[eun E ο, εαν 
BELL DL —— —19 n DN o O 
m Ex a2 YY 42 xy 
o 
O ZN: n 
- 2 DM EN M = 0 (III-3) 
a’ Ze τ xis τ mE 
This expands to yield 
D / M 1 
U - 2U + u y2 ίσου Fly) 
Pl. > ] 9 + [1 = — ip REN uL so 
E XX 2? yy a? Xy 
que rum 
ο ο. -Ἡ : E 
Eo οἱ xt ye " -Σ δεε Y De 
a a a 
Equation (III-4) may be further simplified as shown by 
Landahl [3]. All non-linear terms except the CMM product 
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term can be neglected yielding the following transonic 


small disturbance equation 


Zu 2 um 
[(M’-1) + (y+l)M U. 3 xx 


2M 1 
Ee iuo s Ea cd 
YY > XE as ee 


where a, is the velocity of sound in the free-stream, and 


MES the ratio of specific heats, and Mi=\Mach number. 


B. BOUNDARY CONDITION 

The tangential flow condition requires that the flow be 
tangent to the airfoil surface at each instant of time. This 
means that no fluid may flow through the surface of the air- 


foil and is expressed by the condition 
= 0 On a E) (III-6) 


where 


G(x,y,t) describes the surface of the body 
as a function of time. 


For a thin airfoil restricted to small oscillations, this 


may be written as 


G = y - H(x,t) (III-7) 
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where 


H(x,t) is the function ο σαν πα ENE 
Posicion ος πε απ ο. 


HN, tE) can be προ τοπ for harmonic pitch oscillations as 


Ve 


= d 
Hix,t) = R.P. [a (xx) e ] (III-8) 


where the time-varying angle of attack g(t) is given by 


a(t) = με. --- 
and a, = maximum amplitude of pitcéh oscillation 
Xo is the pitch axis 
ν is the angular frequency of oscillation 
MEC 1 
Pees. = real part of" 


Inserting (III-8) into the flow tangency condition 
auN-6) gives, after linearization, 


$4,040) = v(x,0) = a [U * iv(x - xy] eivt (III-9) 


on y= 0 


This is a condition for the normal velocity to be prescribed 


at the airfoil's mean position y = 0. 


C. NONDIMENSTONALI ZATION 

The terms in equations (III-5) and (III-9) are dimensional. 
For the following calculations it is convenient to use non- 
dimensional quantities. Define non-dimensional time and 


length to be 
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| 


K | 


Aal 


where 


U 


alx 


ajk 


EU 


unitorm velocity pou ποσο νε 


reference length (blade semichord). 


TIT) 


The velocity potential in equation (III-5) may be non- 


dimensionalized as follows. 


Hence: 


e| 
Ι 


Let 


πι 


and similarly for the other derivatives in (III-5), yielding 


2 





2 πω qu cag es Drs m 
El) + (ytl)M ος ie t Ma + MO um 8 


5-12) 
This equation is non-dimensional. 
Pace boundary condition given in equation (III-9) may 
be non-dimensionalized in a similar fashion 
a a ae _ 
9,1%,0) U ne ους (III-9) 
Thus 
Loc 
= E U _ a IVES 
U¿Só,, = EU + ik— (ex - cx] = O (στι ο] 
where 
k = τ΄ - Strouhal number or reduced 
O frequency 
ro ae e E ΠΕ . 
USC yc - NR + ik(x ΧΟ) -= (TII>=14) 
Thus 
D. v(x = l + ik(x - x_)] Qikt (III-15) 
97 = v(x,0) = al ik (x Xo 
Because the final operations are linear in D set a, = L; 
yielding 
B m = = ikt 
ᾧ- - νίκ,ο) = [1 + ik(x - x,)] ε (III-16) 
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The overbars denoting nondimensional quantities will be 
dropped from the remainder of the paper. All further 
quantities shall be assumed appropriately non-dimensional. 
This yields the following final equations 


p 2^ E p 2 _ 
SII o Ot 2M OU M 0 


YY 
ΕΤ] 


and 


ο ο) SA ΙΤ; ΤΠ.“ πο) ] ast (111-168) 


where 


all quantities are nondimensional and 


α = 1 
O 


D. HARMONIC OSCILLATIONS 
In the case of harmonic oscillations, equation (III-17) 


may be simplified still further. 


Let 
Be b + IS 
where 
$ = non-dimensional steady flow potential 
U = non-dimensional oscillatory flow potential 
Bere "real part of” 
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Equation (III-17) then becomes 


- Br h 
(AMÓ Y gg y Μα) νι A O ὁ 


k as 2iMk^y  - ο τ 159) 


For M close to 1, this is a nonlinear mixed elliptic- 
hyperbolic partial differential equation with variable 
coefficients, the exact type depending on ΤΝ and be 
However, because flutter analysis 1s primarily concerned 
with the stability of small perturbations about a steady 
flow, the oscillatory component may be assumed small 


compared to the steady flow potential and therefore the 


product term Uy Vea may be neglected, yielding, 


^^ 


- Ντο) Ib u + ov LU, 


D 
Eos ) Ux ¡es XX’X 


yY 


= 2iM^ky, * M^k^y = 0 (III-20) 
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IV.  LINEARIZATION OF THE GOVERNING EQUATION 


Tne basic flutter equation, ELE ος still a non- 


Mesa, mixed elliptic-hyperbolic partial differential 
equation with variable coefficients and difficult to solve. 


It may yet be further simplified. 


E BASIC SOLUTION 


For M = 1, equation (III-20) becomes 


- - i 2 _ z 
Pyy = τ. Es $i Ud - 2iky, + kp = 0 (IV-1) 


Now assume 


$ ~ Ww = constant (v2 


throughout the interblade channel. Setting 
N 2 


$, Cr*1) = w(ytl) = m (IV-3) 


yields 


nip, 7 yy * ZAR - k^y = 0 (IV-4) 


ο. 





The solution to this equation is found in Garrick and 


Rubinow [5] 


x-my 
í T- l T 2 114 (s-x) 
w(x,y) = v v (s)J, [o Y (x-s) *- (my) ^e ds 
^g (IV-5a) 
LO Ὁ 
and 
x+my 
1 2 i e: 
Ww (x,y) = m / Sn e cm (my) ^]e*^ (5 X) ax 
2 
55; 
LOL EU 
where 
E 9 
MIS Lim «o πε 
Υ-ο Y 
v(x) may be obtained directly from the tangential flow 
Boundary condition, and 
M 
4 
m 
k > ek 
= = lm SS (where this paper 
m^ m employs the 
approximation used by 
Gorelov [4]) 
x, = blade łeading edge, 
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Gorelov [4], has proposed a further simplification. 


Set 


Z = my (IV-6a) 
LAX 


Y(x,Z2) = w(x,y) e : (IV-6b) 


Equation (IV-4) then becomes 


A O 
Xx ZZ 
Bath solution 
X-Z 
Y(X,z) = - = / v(s)3 lu Leer ds (IV—7a) 
X, | 
z 24 
X-Z i 
làs 
τ) IV-7b 
Ez) = E v(s)J. [vo (ese elis zm ) 
m O 
ae 
ZEE 
where 
Ed A .. 
iS Ae lim 2,00, 2) 
σου 


v(x) is obtained from the tangential flow boundary condition. 
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For a thin body immersed in the flow, the solutions 
moe Y - 0, Z > 0, and y < 0, z< Q apply above the body 
along left-running Mach lines, or below along right running 


Mach lines respectively. 


B. BOUNDARY CONDITIONS 
1. Flow Tangency Condition 
The boundary condition comes from the tangential 
Bow condition, (III-18) 


v (x) = [l + ik(x - κ.) | (IV-8) 


2. Upstream @onaltion 


The final linearlized equation is a hyperbolic differ- 


ential equation with boundary condition 


vix,y) 0 (IV-9) 
when 
X=“ X, < | my | 


for the solution shown in equations (IV-5) or 


πα“, τι 
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V. PROBLEM FORMULATION 


A. CO-ORDINATE SYSTEM 





Assume the geometry shown above. Both blades are thin 
alrfoils of semichord c. All measurements are non-dimensional, 
normalized to c. The (x,y) co-ordinate system has its 
origin at the center of the reference (lower) blade. The 
(x, Yi) System is centered at the middle of the adjacent 
Mapper) blade. The origin of the (x, Y}? system is located 
at (0,p) in the reference system. Generalizing this conven- 
tion, the same symbols shall be used for the same quantities 
on both blades. Where discrimination is required, the quantity 


associated with the adjacent blade will be marked with 
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superscript i the quantity associated with the reference 
blade will be either unsuperscripted or marked with a superscript = 
Each blade is assumed to perform a small amplitude harmonic 
oscillation about its mid-chord point. Both blades are assumed to have 
identical reduced frequencies, k, and the motion.of the. adjacent 
blade lags that of the reference blade by a phase angle c. 
The blades are immersed in a uniform flow from the left 
at M= 1. The objective is to determine the oscillatory 
pressure distributions and aerodynamic forces generated by 


the blades' oscillations. Cascade solidity, > 2/p.- 


B. BOUNDARY CONDITIONS 
l. Upstream Condition 


y 20 whenever 


Ket ΠῚ 


(V=) 
and 
x + 1 < my, |; simultaneously 
2. Flow Tangency Condition 
Along the reference blade 
lim y, (x,y) ^ (1 + ikx) (V-2a) 


Along the adjacent blade 


87 





1σ 
e 


lim us (X, y) = (1+ ikx) (V-2b) 
m 


y, >0 
where 


σ is the phase angle between the blades oscillations 


fee BASIC SOLUTION TECHNIQUE 


Assume that the unsteady potential, y, may be written 


as the sum of four components 
ψίκ,γ) = ὁ ίκ,ν) + P Guy) * e Guy + 0 (x,y,) (V-3) 


where: 


$ = potential due to the reference blade 
alone, known from equation  (IV-7) 


$ = potential due to the adjacent blade 
alone, known from equation (IV-7) 

y? = interference potential required to satisfy ~ 
tangential flow condition along 
reference blade, unknown 

Ta = interference potential required to satisfy 


tangential flow condition along 
adjacent blade, unknown. 


This total potential must satisfy the tangential flow 
condition at the plane of both the reference and adjacent 


blades. Thus 


O - T A O B JE = 5 
5Ο + τ EN y + a! p) (V-4a) 


o JE x) 
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at the reference blade, and 


O = ÉE 5 Ο x 1: . 
6 (x,y=p) + Bog ze ee E el 


= (1 + ikx) e'? 


at the adjacent blade. 
But from the unsteady potential solution for a single 


oscillating blade one has 


97 (x,y=0) - 1 * ikx 
and 
ik _ . ic 
y A Y A e 
E E 
Thus 
= ( EU yt (X =-p) + Vo (x y=0) = 0 
E ds πα ya 
along the reference blade, and 
6° (x,y=p) + vo (x,y=p) + UL (x,y,=0) = ο 
ed Y Nr l 


along the adjacent blade. 


From equation  (IV-7) 
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(V-4b) 


(V-5a) 


ο 


(V-6a) 


(V-6b) 





x-my 


O JE j = 
ὁ (κ) τα [| —vÜ(y pot(x-s)?- (ny) 21 ME as 
πι 
σι 0 (V-7a) 
x+my 
il = 
= = f vy pulis) - iy) ^16 ^ 579) as 
zZ 
y < 0 (V-7b) 
where 
v? (s) = 1+ iks 
= k/m 
m = (y+l)w 
| _ κ΄ (1-π΄) 
τ à 
m m~s 
W = mean value of oy in the channel 
X-my4 
IE EET l : 2: DN Ss 
6 (x,y,) = = Í a aa “Je ds 
' 
Μ 20 (IV-8a) 
x*my, 
PO i a^ 2. _ iA (s-x) 
E I v (s)J, luV(x-s) - (my) le ds 


EL 


Yi < 0 (IV-8b) 
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where 

Y ο - (l1 + iks)el? 
Henceforth attention will be restricted to the flow within 
po channel, 0 <y <p , -p < Y< 9 leaving 3. anda 
(IV-8b) as the governing equations of interest. 


The two interference potentials are assumed to have 


forms identical to (IV-7a) and (IV-8b). 


Set 
x-my 
ον αλ {ο 
Vv) -- = u^ (s) J, [à (x-s) ^- (my) 167 “πολης 
-1 
y > 0 (V-9a) 
Schm 
2 2, ià(s- 
ντ) = = J ul (s)2, Ly (x75) - (my ) 16^ (979 as 
= (V-9b) 
YO 
where 
u? (s) and mS are unknown functions to be 
determined so as to satisfy equations (V-6) 
Substitution of (V -7a), (V -8b) and (V -9) into (V -6) 
yields 


ul) + pl aya) * 0 


E y DIS (us poa 
Yi Y i 


4l 





ut (x) + p, (x, y=P) T 6. (x, y=P) = 0 (V-10b) 


Recalling Gorelov's transformation discussed in 


r 


[4] and shown in equations (IV-6) above, set 








LÀ 
$?(x,z) = 9 x.y) e” 
l zs . pol 
ᾧ (κ, Ζ1) ο φ (X,y4) e 
LÀ 
YOlx,z) = v lx y) e * 
IL 1Ax 
TU(x,24) ψ (x, y) e 
where z = my 
27 = p 
Then 
qui x 
= we) to (X,2477X4) t * (xX,217-7X,) = 0 
= 1 1 
(V-11a) 
and 
LAX 1 O O 
ze x et aX TB y EO 
m Ζ Z 
where 
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To employ the collocation method, assume that ut (x) 


O ! 
and u (X) can be approximated as the sum of a set of elemen- 


Bey functions f, so that 


J 
n 
u^ (x) X ) 9: f j (x) (V-12a) 
I) 
n 
l 
u(x) z J T e bn (V-12b) 
1Ξ0 
where ES TO wen Xx < x,-1 


Note that here both u? and m are expressed in terms of the 
Same elementary functions, fj. 
Because of the slightly supersonic nature of the problem 


observe that u°(x) = 0 and al (x) IS IE 


Equations (V-12) may now be rewritten as 





1 
1AXr,O 9 L ο τοι AS 
e 193. (1) - 3z, f 163£4 (S) Js lu (x-s)'-z1]e 85 
xl 
Xt*24 
m ig. o a OE EAS 
5 Nie f (1+iks)I  [uV (x-s)"-z]]e "ds 
1 
-1 
Bu Lc (NES a) 
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xpa l . ο. ο 2 lAs 
e JOE, (x) το / 18 5£ ;(s)Js lu (x-s)^-z]e"^"ds 
el 
DT 
5 | 
EE i (1+1ks)J [wY (x-s) f = z^] e^ Sas 
O 
SL 
απ νο ο κ. ο πο) 
where 
οι ο του XU T 


This simplifies to 


AS 





AMS + jet ? J GS lu (x-s) ^-221e? ds) 


OZ 
l 


XTZ 
1 
ng 9 ES 
en n (leiks)J, [uN(x-s) -z7]e ' ds 


E 


MEC M (V-14a) 


and 


PAX ca Or ð Lr ee irs, 
e 19.6. (κ) -}9. {πᾳ ji f(s) J, lw (x-s)^-z^]le S 


2 id 
- 2 f (itiks)J_[uV (x-s)?-2*]e**Sds (V-14b) 
92 O 
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aZ = x. where 
eS = 0 for πο τν 


emang the indicated differentiation yields 


A J, DoY (x- s) ^-x2 10x, e^ "Sas 


iÀxvao l 
D P. Bor 
e 195€; (x) + 28, / ues) 


J πρ A: 
Zu (x-s) ^-x* 


* 


= E (xx) e Ts) η 

; um. ο. τω Bea. 

f E AA ee 

p 4 x-s) ^-x2 
ià (x-x.,) 
- [l+ik(x-x,)]e A (V-15a) 
and 
= INS 
i νο J. I0 (κ-θ) ?=x7Jux e+ ds 

ei^Xpyele (κ) 1 Jet f 2S 


J I J πο το 
αμ (x-s) ^-x7 


- ics "d 


ωχ J 


af DEED 
cm "ds λος κ) 


: nr NE 





= -= Í (l+iks) 
Zi E 
(V-15b) 
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Gorelov's formulacion, equations [2.6, 2.7, 2995 and 2.9] 


Seles, can be obtained directly from equations (V-15) by 


substituting 
f; (x) = COS Ja EOS DN 
n 
v? = J) u? cos jn 
15ο η 
n 
y? =) 2 cos 3n 
j=0 
where: 
ee = cos (=x) 
NES - σος Γάτος 


In comparing the two systems care must be taken to note the 
differing symbols and coordinate systems. The corresponding 


quantities are: 


Here απ, 
O a 
M = Yoo” “io 
O 1 
RO "3 9oc' "ic 
n 7 cos t(x) n = cos *(1-x) 
M cos™t(1-x,) Dij os cos | (1-x,) 
2, Zu Y: i 
j σ 
σ Y 
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Eu x«1;in [4] 0 < x « 2. This transformation 


accounts for the differing definitions of n. Making the 
substitutions results in the system 


n 


) (97 [cos HE (17954! tos 31) 
1-0 


De 


+ Ur if a, 


Zool 


1 at _ . πλ κος] 
+ 9; [cos an (1 $:9)cosjn,le ) 


n 


Rex 
E. D 9 Be -*_1XS 
== Lu: / δρ οἱ (x-s) 71] cos Jne ds 
j^0 τι 
-- ul cos nee } (V-16a) 
at 2] = 7X, 
πε’ 
and 
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2087 LA i al 
32, NU (x-s) -2,] [cos jn - (1-8, )cos jn,le 


AS 


ds 





E ; 
L (97 [cos jn - Ue cos jn,] 


ux 
ο 9 NE LE | 
- D / pali xs) 7-2 ] [eosjn- (1-6, )cosjn,le~ “ds 
poe 
O pa Js 
D z zuge : * 
"E ad SS jn,le ) 
n κ. 
_ O 9 «/ D ου x 
= J tu. f SE Jg lo (xs) -2 ]cos jne Sas 
3=0 -1 
η (s σος ye AX o) y welch) 
where: 
n= arc cos (-x) 
io = amemcos (1—x,) 
N = are cos (=s) 
n = arc cos (-xłtx,) 
ib for rj = 0 
oe = Dirac 6 function = 
93 O O 


Given the change in coordinates and notation, this system 


is equivalent to that shown in [4]. 


48 








This was the system programmed for computer solution. 
Because the function, Ez (x) + is unaffected by the differen- 
tiation with respect to y (or z) the exact form used need not be 
specified, so that the system shown in (V-15) may be programmed 
with f undetermined. A subroutine may be written to return 
the function desired and the remaining program left perfectly 
general. In the program developed with this thesis both the 
Gorelov functions shown above and the Legendre polynomials 
were employed. All the integrals may now be evaluated at 


n+l points, X;¡» on both blades in mp-1 < ο l and the 


resulting linear system solved for 93 and 95, 515,25 ο δις 


The interference potentials may be constructed by taking 


x- zZ 
O 2A -1 Ο τα 28 2 DAS 
ee) = — Í $ oE (x) 1J, fu (x-s)^-z^]e '"ds 
z- 20 (V-17a) 
APS, 
l iL l ο. ELAS 
y (x, Z) = / [203 £ G0 12, Do Y (x78) -z,le ds 
-1 
Z4 < 0 (V-17b) 
where 
£.(x) = 0, for aki xul 


Once the potentials have been calculated as outlined 


above, the surface pressure may be calculated using the 


relationship 





C, - -2(0, * ik) (V-18a) 
ER Πες. (V-18b) 


Because all the plates are assumed to be in steady oscillation 
with uniform phase shift, o, between neighboring plates, 


then 


vt (x) - v? (x) e?? ' ut (x) = hero eas meee y) Sexy) 


in this case 


| 
ο τα | [v (x,*0) + NU (V-19) 
Eh 
where 
X 
¥ (xX t0) = = J [(L+iks) &u? (s) 12, (o (x-5) 1e] ^ (9) as 
E 
ve X X ------------- 
τ q Í | (L+iks) +u? (s) 1J, [uV (x-8) 2- 0 21e ^^ as 
=l 
where x, = Mp. (V-20) 


Results from this approach, in the form of values of Cta 


for k=0.1, at various values of w are presented in the results for 
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approximations based both on Gorelov's formulation, and on 
the Legendre polynomials. 
D. COLLOCATION SOLUTION OF THE POTENTIAL EQUATION EXPANDED 

FOR SMALL k, 

In order to provide a partially independent check of the 
results of the main program, the Gorelov function represen- 
tation of the collocation solution was expanded for small 
Miñana solved at two collocation points, n= 2. The resulting 
potentials, and partial derivatives with respect to x and 
y were then used to replace the corresponding numerical 
routines in the main program. The output resulting from the 
approximations were compared with the purely numerical results 
obtained from the computer program. 

1. Solution For The Unknown Potential Coefficients 

The basic system of linear equations used to determine 


the unknown coefficients is 


MA O l JE E - rm eu 
m uer UM Ἔ z^ = 0, 2 0, 2, S mp 
l 
Ru) 
- iu ο ον ο, 2 5 09 uz EE EID (V=21D) 
m 2 Ζ L 
where 
0 = uP(x) = ul(x) when x « -l*x, 


ο] 





otherwise 


n 
OW. O : . ο 
mx ἡ 9,(cos jn - cos jn,) + 69 
ge 
n 
JE E 
u Ξ ) 95 (cos Din ces Tr a? 
J=1 
where 
n = arc cos (-x) 
are cos Ml==x,) 
Thus, for n = 2, the system becomes 
eL ^X(g9 + 97 (cos ου. οσα ης) } 


KFZ ; 
2 2 M MEE c 
97 O A 
κ. 
XT*Zi | 
ον f ulls)s_luVix-s)°-z’le as = 0 
95 O Hr 
1 
a 
Z1 = “Xy (V-22a) 


along the reference blade and 


-- 





NZ 


acl E _ m9 ο ec, ts 
EN 2, (cos n = cos n,)} - 5, I (s)J, Lo (x s)^-z^]e "ds 


ll 
pm 
LÀ 
-Š =f δώσ [ωγίκ-ο))- αἼο) δας. (Υ-225) 
Kl z-x, 
along the adjacent blade 
where α- (5) = = + 97 (cos me cosin) 
ΙΙ» Ξ gt + 61 (cos n = cos n4) 
O T 
v (S) = 1+ iks 
1 u πο 
v (s) = (1 + 1ks)e 


For k sufficiently small, this system may be further 


simplified by the following approximations 


J [w "ΟΜ. (Y-23a) 


è? 


5-2... o a (V-23b) 


ο 
RR oe wee προ. ee ο» πο αλ 
j 2 8 
115 _ l + iàs = 0 (A*s Je alan = 
where 0 (u^) means "of the order of magnitude of ο. 
ο x «Ll, ΛΔ S S xXx. 
2 2 
ne eae ο; _ K D ) 
m 


DS 





The interference source distributions may be replaced by 


O - O ο. . 
u (s) = E ΞΕ 94 ( SEX cu 
ih u l 08 x 
ο ο = gs + 94 ( So mu) 


If higher order terms are neglected, the result is a system 
linear in k and λ 


KEZ 





O ee ; 
(l*iAx) [02 E 4 í (1+iks) (l+iìs)ds 
Dd n 
XZ 
3 i l - 
+ — / [97 + 07 (-s-1+x,)](1l+iAs)ds = 0 
921 O l 
xvm] 
za + (V-24a) 
KFA 
EMEN sbox-12x0] --— f  (1-1ks) (l+ids)ds 
O H = 92 
-1 
xm 
Q : E 
E = f 182 + 0, (-s-1+x,)] (l+i1s)ds = 0 
QT 
. (V-24b) 
Ze 
2 
Product terms containing (kA) = — may be neglected as 

m 


of higher order in k, yielding 
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i O O 9 
(1*iAx) [85 + 8) (-x-1+x,)] + 5 f [1 +4 (Atk) s]ds 


X+Z 
9 1 i 
* 9z. f t [87 + 8) (-s-1+x,)] 
pe s 
+ iàs[o} + el(-s-lex,)l)ds = 0 
παν»: (V-25a) 
22 
(1+idx) [82 + 87 (-x-1*x,)] - fJ  U*iQsk)slás 
EXE 
κ-Ζ 
9 O ο ο ο 
En f MA Es iAs[9, + 98, (-s-1*x,)l)ds = 0 
pel 
Zo cL 555) 


Evaluating the indicated derivatives yields 


(1*iAx) [00 T Ur T ME VE E iene 


+ tol + el (2x,=x-1)]+12 (0x4) [97+0] (2x,2x-1)]) =0 


(V-26a) 
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(1+idx) [05+01 (-x-1+x,)] + [1+i (Atk) (x-x,) ] 


η i £) (x-x,) [8048] (2x,-x-1)]} = 0 


(285) 


TRUS: 


9^ in) E NL E + ο ο ο νο) 


* ol(l id(x-x,)] + 9l[(2x,-x-1) * iA (x-x,) (2x,-x-1)] 
E Bo Lo cx] QUEE NEW 
Ti EEUU) + olr = = z 
5 ΙΙ 


E το ος κλίκ.) (2x, ΚΤΠ 
ο ο ο ο οκ ο ο. Το) 


This system may be solved at two points, x, and X» ; rer 


O O JE 1 
E θη, or and 97: 
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2. Calculation Of The Potential 


The potential is given by 
W(x,y) = "(x,z)e b 


where 


x 


ιο) = = 


E 


= 


/ [v^ (s) *u? (s) 13, [o E τ; 


(V=28) 


i = J [v3 (s) u^ (5) 12, [o p 


u-(s) = u'(s) = 0 που alle: SM 


Thus 


X 
BExZ) = - z f vo (s)3_ [o (x-s) le Sas 


zi 


x 
f. u? (5)2, [oy (x73) *e 


poc 


LAS aa 


Zn 


XU. 


+ Es / v3 (s)J, [u 7 


ds 


* lÀS 


+ = / ut (s)J, [u (as) ee 
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(V=29) 





Making the same small frequency approximations as 


in the previous section yields 


X 
¥(x,2) = -= f vs) (1+ids)ds 
-1 
x 
- = f u“ (s) (1+ils)ds 
ore 
1 Xe 
dE : 
.- 1 vo cnc sss 
-1 
Kein 
1 i ; 
+ = f Un (Sues yids (V-30) 
el 


From the general formulation 


Yy= 6 +6 +V7 + y Vz) 


Thus, along the reference blade 


X x 
-md-(x,z=0) = f v°(s)(1+ids)ds = fJ (l+iks) (l+iAs)ds 
ET =| 
s? x 
= f [1+i(k+1)s]ds = [sti (k+A) =] (V-32) 
E zu 
= + OH ic iE 
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O 
$ (x,z=0) = - ={ (1+x) + i (Em ESSE (V-33) 


By inspection 


O 





1 
9 Guz,-x,) = Sf (lex-x,) * i825) [(ex,)2-1])  (v-34) 
X 
O O . 
=m¥ (x,z=0) = f u~(s) (1+ids)ds (V-35) 
pe-l 
X 
- f [95 + ef (-s+x,-1)] (1+ids)ds 
τα 
Χ 
E f 95 (1*iAs) + 8] (-s+x,-1) (l+iAs)ds 
xm 


- 921 (alo) 5 x^ -x 282x471) ] 


X 
2 f e 


1 


O 


O [xx 1412 (só+sx,=8) ]ds 


= 89 [Gclox, H5 t2, 1) 


O s2 gs? s s? . 
u -— 1 - 
sa 93 t MSX, S)+iA ( uta cse] ana 
* 
(V-36) 
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u O A LA 
= 9,l(xt1-x,)+5-(x 7x, *2x4,-71)] 


(x,-1) 4 
1 2 * 2 ΓΙ ΤΙ 


3 x^x, ο. οὐ 


| E En ΠΡ 
E. uu oen el 


Ξ 99 [ Gcel-x,) 5 a x 12x, 71) ] 


3 
, x 2 (x,-1) X 2 1 
+ id[- Gea mgr zi Se 
3 2 
Xx 2 Xx Xx 1 
νο το. τοι; = X le 
= rea) + 2x3 zx, - 001 
2 2 
ο, x 5 
- 971 ud x 1) 2 
3 2 3 2 
EU. oux Cum ee ids 
ο mm 2 mu 2. * gli 
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e .. MO πλ, 2 2 
M = > o + — x ου... ΠΠ] 





2 2 
+ 62(- 5- e x(x,-1) - MD 
DE PR e 
* X yl) 
+ iAl- — + 7 Eo —]1] [νη 


JE 
my (X,2477X,) NE ccevalmumsed by Substituting x=x, for x τη 


the expression for πι” and exchanging θ- and A τος o^ and 
ο 
a 
my (x,z,7-x,) 
LÀ 
97 CL (rox, ) 17x, + SAU (x-x,) * = x^ B AE ΙΤ 
1,- (x=x,)? E 
πα ο. 
mesi (x-x,)> B x? (x4-1) 2 c (V-38) 
3 2 6 
= a Educ )ct SAL (x°+2xx,+%5) = x? E OE 
2 2 
lord et e T 23 d , 
E 971 Ds, (x πι. l) 
x? 2 2 κὶ 2 2, (xe 
+ id[- F E ee t xs + (x LR) πο 
(x Ls 
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4 


» al πλ AO 
= ο ος) ο ο ο ο εν οκ ο UN] 
ol; x x? T 
+ n > e ο > + x(x-1) + >(1-3x,) (x,-1) 


B : 2 2 
EEUU x m κοκ; T e 


Bs 
+ Axe ee l 
2 


y? girs 


i 
(82 { (x+1-2x,) m x? A ee] 


| 
ES 


E θτ{[- = ax = em xx? PAZ I 5] 
O Ce) fed x al 
Em ca plas ste om 
(V-39) 
A comparison of the results for the full program and the 
approximation is given below for k = 0.01, w= 0.05, 0 =P, 
NEN yielding »- = 6.1x10 7, X s 0.083 
Bull program Approx 
x = , 1285 
È ο -3271 25.379, 5491 
Oa =8. 634, „3011 -8. 66, 73791 
xX = .5643 
È Eom, .8091 -9.804, 1.000i 
Za, 6921 212.2724 :7608i 
Χ 
Co +31.658,-3.1032i τυ alo Sl 
a α 
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VI. RESULTS 


The collocation method was used to solve the partial 
differential equation resulting from the Gorelov approximation 
of transonic potential flow in an unstaggered cascade. The system 
was solved using both the spanning functions proposed by Gorelov in 
[4], resulting in the equations (V-16); and the. Legendre poly- 
nomials, resulting in equations (V-15) with f. replaced by 
the Legendre polynomial, ZH The resulting values of er 


for k = .1l, T = 1, o = 1 and seven collocation points on 
each blade are presented in figures VI-2, VI-3, and VI-4. 
Figure VI -1 presents a diagram which is useful in 
commenting on the other results. This shows the location of 
the collocation points and first three interference reflec- 
tions as a function of w expressed as a percentage of that 
portion of the chord subject to reflection. The collocation 
points are equally spaced throughout this interval, 12.5$ 
from the leading edge of the interference zone, 12.5$ 
between each pair of points and 12.5$ from the blade trailing 
edge. The independent variable, w, is plotted vertically 
5ο that the dependent variable, percent of chord subject 
to interference, may be more conveniently visualized along 
the blade. (The curves are not precisely linear, but are 
very nearly so in the range shown.) 
Figure VI -2 shows the Co calculated with k = 0.0 in 


comparison with the results obtained from Ackeret theory. 
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Agreement is good where there is no reflection and the 
portion of the blade subject to interference is affected 

by a constant interference potential, w ο oorheewuout 
the rest of the curve the results calculated here oscillate 
above and below the theoretical values. This appears to be 
due to the discrete nature of the approximation used in the 
collocation method.  Rarely is the fraction of the chord 
subject to interference reflection equal to the fraction of 
the collocation points which feel it. Where the collocation 
point fraction lags, as near w = 0.6, the collocation results 
are lower than those due to Ackeret theory. When the 
Euocation point fraction leads, as it does for w < .04 

and briefly for w = .08, the the collocation results are 
higher than those due to Ackeret theory. The fault appears 
to be an intrinsic feature of the small number of points 
Sampled. This results in a set of coefficients similar to 
those which would be obtained from a generalized Fourier 
series based on the integration of the Taylor series 
expansion about each point. This obviously cannot be a good 
approximation when both the function and its derivative are 
discontinuous at the reflections. 

Figures VI -3 and VI -4 show the results of using 
Legendre polynomials and Gorelov's functions as spanning 
functions. The results for both formulations are identical. 
Gorelov's results are presented for comparison. Agreement 


is good for w > 0.05 except for an anomalous point, marked A. 
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It is believed that this anomaly is due to the location of 
Bue first reflection just ahead of the last collocation 
point (cf. "A" on Figure VI -1). This will yield a very 
small contribution from the reflection potential to the 
linear system from which the collocation points are 
determined. The resulting system will have a large dynamic 


range and may be ill-conditioned. 


The discrepancy between these results, and those in 
[3] for w « 0.5 is still unexplained, as is the outlying 
τις for w = 0.5. 

The discontinuities in the imaginary results are believed 
NO Be due primarily to the reflection/collocation inter- 


action explained above. 
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VII . RECOMMENDATIONS 


There are two recommendations to be made about the 
techniques used in the collocation method, and a new area 
in which it might be employed. 

The program developed in the course of writing this 
thesis employs adaptive Simpson's integration to calculate 
the elements of a completely determined system. This system 
is solved to provide the coefficients of the spanning functions. 
Two improvements may be made: 

l. The Simpson's integration scheme may be replaced by 
a Gaussian integrator. Experience has shown that several 
thousand function evaluations are required by the Simpson's 
integration routine when Ca, is to be evaluated for small 
w. This entails large amounts of computer time and leads 
to increased accumulations of numerical error. Use of Gaussian 
integration would probably improve both of these characteristics 
with little loss of accuracy. 

2. The present program treats a completely determined sys- 
tem of dimension 2n+l by 2n+l, and then solves that system to 
find the collocation coefficients. This procedure has worked 
satisfactorily in this thesis, but may not work as well at 
higher frequencies where the final linear system of equations 
may’ be ill-conditioned. As an alternative, it is recommended 
that the boundary conditions be applied at more than n points, 
say m points, where m is twice or three times as many points, 


and that the least squares technique be used to determine the 
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the collocation coefficients which give the minimum square 
error over-all. This may be thought of as "sampling more 
data" in order to get more information about the unknown 
function. The present program could be easily modified in 


this regard by replacing the spanning function matrix, QlZINT, 


by a new matrix of the form 


QlZINT' = xix 


where X is the new m by n+l (m > n+1) matrix, and replacing 


the present right-hand-side vector, Ql1COF with 


olcoF' = xy 


where Y is the new m by 1 right-hand-side vector. An 
alternative would be to employ a prepackaged statistical 
linear regression routine after either modifying the routine 
to accept complex data, Or transforming the present system 
into a larger system of real numbers only. 

The new area in which the collocation method might be 
employed is the calculation of the potential flow about a 
staggered cascade. The method coud be employed to calculate 
both the potential in the channel and above the upper blade. 


The program presented has been designed to enable the 
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calculation of flow within the channel of a staggered 
cascade. Unfortunately, there was not enough time to extend 


the study to this case. 
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APPENDIX A 


PROGRAM DESCRIPTION 


This section describes the computer program used to 
calculate the interference solution to the Gorelov lineariza- 
tion for unsteady transonic flow in a channel. The program 
written in IBM Fortran IV with the basic structure outlined 
Mis tcevens [5]. The basic points are: 

l. Organization of the program into small subroutines, 
each of which performs a specific task. 

2. Transmission data to and from subroutines via a 
formal parameter argument list. No common statements 
are used. 

The end objective is code which is both easy to modify and 
maintain. 

Each subroutine is designed with optional diagnostic 
Beamneing of its input and output. This is controlled by the 
Beapamecer IPT. The diagnostic output is printed (only) if 
IPT > 0. Each routine accepts IPT, sets IOT = IPT - l, and 
then passes IOT as the print parameter to routines it calls. 
By this method, diagnostic output can be "cascaded" to any 
desired level. Large initial values of IPT should be avoided 
because of the spectacular amount of output which can be 
generated by the double integrals within QIDCOF. 

1. Main Program; including subroutines READ and ABSA. 

The basic structure is given in Figure A-1. MAIN 


calls READ to read input data and then ABSA to calculate the 


p 





MAIN 


ie nn 
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QlZINT 
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Ως ΡΕ 
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OTERTE 





Figure A-l. Program Hierarchy 
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collocation points. The version of ABSA shown evenly spaces 
the collocation points across that portion of the blade 
subject to interference.  ABSA may be easily replaced if 
different point spacing is desired, or if additional points 
are to be added for an overdetermined system and least 
Squares approximation. 

2. QSONE This subroutine controls the actual potential 
calculation. It performs no calculation itself, but calls 
subordinate subroutines where the calculations are actually 
performed. The calling hierarchy is shown in Figure A-1. 

3. QLZINT This subroutine calculates the linear equa- 
tion system arising from the boundary conditions. Hierarchy 
WS Shown in Figure A-2. 

τους matrix output is carried through QlINT.  OLZINT 


calls the following subprograms 


a. QLICRBP returns the value of 6, 
b. QICABP returns the value of oa 
C.  QICRIP returns the values of ve 
Χ- 2 
9 | B 2-2 
= / £5(s)Jg IuV (x-s) -z ]ds 


-1+x, 


where fi is one of the set of i elementary 


functions, j=l,n 


1 
d. JQICAIP returns the values of Y, 
1 
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XT 


Z 
a J ee as 
BUOOPESEU Ex 


where = is one of the set of elementary 
functions. OFFSET is a.parameter 
included to facilitate program 
conversion to a staggered cascade 
e. PLGNDR 1s the subprogram which returns Se 
the elementary spanning function. No other routine contains 
ENpircit reference to the spanning function. This facili- 
tates easy replacement of the spanning functions should this 
be desired. | 
OICRBP and QICABP in turn call QlCWXM and QlCZXM. 
QLCWXM computes end-points and then calls Ql1CZXM, a complex 
integration routine based on SIMP by Shampine and Allen [6]. 
QICAIP and QICRIP call QICWSM and Q1CZSM to: perform the 
integration.  QlZINT passes the constant matrix to QlCOEF 
in the array QlINT and the right-hand-side vector in the 
παν OICOF. 

4. Q1COEF This subroutine employs the IMSL routine 
LEQTlC to solve the linear system received from Ol1ZINT. The 
resulting coefficients are QlABCF for the adjacent blade and 
QIRBCF for the reference blade.  LEQT2C, the high precision 
complex IMSL routine may be directly substituted for LEQTIC. 
OLCOEF may be rewritten to employ the generalized inverse 


required for least squares approximation 
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Ah = 
ΚΠ) x y where x 


_ QILINT 


QUEDE 


KK 
I 


after first performing the multiplication necessary to replace 
ΙΤ and QICOF with the proper matrix products in the call 
NO LEOTIC. 

5. ΟΙΡΡΟΤ This subroutine calculates the potential, 
$, and $,: at each collocation point along the reference 
blade but only if QIPPOT receives a value of IPT > 0, re- 
E "Hj Pr > 2 on input to the main program. I£ IPT < 0, 
then the subroutine is exited before any calculations are 
performed. This subroutine is most useful for debugging 
ΠΡ and QICOEF.  QIPPOT calls QOlPABP, QlPRBP, OQIPRIP, 
ENDEIP,OIDABP, QIDRBP, QlIDRIP, and QIDAIP, all of which will 
be described in the next section. 

6. QIDCOF This subroutine calculates the dimension- 


Mess coefficients of lift and moment; Co.r Cp. Its 


Q Q 


internal hierarchy is shown in Figure A-3. 

a. Ql1DCL calculates the nondimensional complex 
coefficient Cy. 

b. QlDCMcalculates the nondimensional complex 
coefficient Cy,» 

C. QIPRBP and QIPABP calculate the potentials due 
to the reference and adjacent blades respectively. 
O1PWXM and Q1PZXM are called to perform the 


actual integration. 
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QIPRIP and QlPAIP calculate the interference 
potentials along the reference and adjacent 
blades.  QOlIPZSM is called to perform the 
integration. 

QIDRBP, QIDABP, QIDRIP, and QIDAIP correspond 
exactly to subroutines above except that the 
value returned is the partial derivative of the 
potential with respect to X. QI1DWXM, Q1DZSM, 


and OIDXSM perform the co-reponding integrals. 


6. Program Listing The program listing shown below 


incorporates the Legendre functions as spanning functions. 


Listings for a subroutine employing Gorelov's spanning 


function and the linear approximation program follow. 
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8. Linear Expansion Program 
The program based on the linear expansion for 


small k is shown below. 
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